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Abstract 

- Context Mechanical wood properties are increasingly rele- 
vant for structural applications and are influenced by growing 
space availability. Douglas-fir (Pseudotsuga menziesii (Mirb.) 
Franco) has an increasing market share in Europe and is 
mainly processed to sawn timber. 

- Aim A sample of 164 thinning trees was taken from two 
Douglas-fir long-term forestry research plots in Germany. The 
end-use quality of about 2,000 side and center boards was 
analyzed as a function of initial plant density (1,000, 2,000, 
and 4,000 trees per hectare) and log position within the stem. 
¢ Methods Sawn timber quality was described by knottiness, 
density, modulus of elasticity, and strength. Explanatory pa- 
rameters were radial position, longitudinal position, and initial 
plant density. All boards were strength graded visually and by 
the grading machine GoldenEye-706 using both X-rays for 
detecting densities and size as well as position of knots and 
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laser interferometry for detecting eigenfrequency (DIN 4074, 
DIN 2012; EN 14081-2, CEN 2010). 

- Results High plant density led to better mechanical sawn 
timber quality. Significant differences were especially ob- 
served between 1,000 and 2,000 trees per hectare. The yield 
of machine strength-graded center boards of strength class 
C24 increased from 50 to 89 % at low and high initial plant 
density, respectively. 

¢ Conclusion Foresters are able to improve end-product qual- 
ity by controlling planting density in particular. The round- 
wood price that foresters get should be based on the propor- 
tion of higher strength classes within logs to give incentives 
for a more quality-oriented forest management. 


Keywords Wood quality - Knottiness -Dynamic modulus of 
elasticity - Strength grading - Tensile strength - Board position 


1 Introduction 


Forest management of softwood has moved towards more 
integrated approaches that consider risks such as snow, storm, 
and bark beetles (Griess and Knoke 2011) and imply wider 
initial spacing (Smith and Reukema 1986; Hein et al. 2007). 
Silvicultural guidance to the forester usually aims to maximize 
volume yield and takes into account the multifunctionality of 
forests but pays little attention to the actual properties of the 
material produced. There are concerns that the relevance of 
timber quality might deteriorate in the future, and, as a con- 
sequence, a better understanding of the processes that control 
wood quality is required. 

The second National Forest Inventory (NFI2 2005) has 
quantified the current area of Douglas-fir forests in Germany 
at about 180,000 ha, equalling a portion of 1.7 % of forest area 
in Germany. The share of intensively managed forest stands of 
Douglas-fir has increased in Germany within the recent 
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decades, particularly in the 1960s and 1970s. Due to this 
significant increase in planting, currently, many trees are 
considered as tree age class 2 (21-40 years) and class 3 (41— 
60 years). As thinning takes place at these stands, an increas- 
ing amount of logs will come into the market. Most wood will 
be used as a structural material. The assessment of Douglas-fir 
timber quality in relation to forestry measures is increasingly 
considered as important (Hapla 1980; Sauter 1992). Naturally, 
the variance of the quality of Douglas-fir (sawn) timber is 
large: Douglas-fir coming from the original old-growth forests 
of western USA with low annual increment is superior with 
regard to wood quality compared to the fast-growing trees of 
modern-managed short-rotation forests (Constantino and 
Haley 1988; Barrett and Kellogg 1991). There are indications 
that the available forest resource of Douglas-fir timber has 
changed and has a negative impact on its strength and stiffness 
values (Dahlen et al. 2012). 

Within the upstream section of the production chain, forest 
managers may take fundamental decisions for wood quality of 
Douglas-fir that are related to (1) site conditions (Kantavichai 
et al. 2010); (2) genetics, breeding, and selection (Vikram 
et al. 2011); and (3) silvicultural treatment (Robbins 2000; 
Weiskittel et al. 2006; Hein et al. 2008; Newton et al. 2012). 
Plant density in particular is a cost factor and has a significant 
influence on wood properties. For other softwoods, plant 
density determines the management of optimum wood quality 
at a considerable extent (Brazier and Mobbs 1993; Robbins 
2000; Kang et al. 2004; Fujimoto and Koga 2009). 

Wood processing may control the physical and mechanical 
properties of sawn timber through board position and sawing 
pattern (Lindstrom et al. 2009; Ovrum et al. 2011). Besides, 
timber quality of Douglas-fir can be nondestructively predict- 
ed by acoustic tools at the level of standing trees (Wagner et al. 
2003; Amishev and Murphy 2008a), of logs (Amishev and 
Murphy 2008b; Rais et al. 2014), and of sawn timber (Dahlen 
et al. 2012). For structural use, timber quality is crucial and 
usually defined using indirect indicators such as a quality 
grade that mainly depends on knottiness or direct indicators 
such as density and stiffness (modulus of elasticity) which are 
predictors for strength. In particular, stiffness is a very impor- 
tant indicator for log and sawn wood quality for Douglas-fir 
(Dahlen et al. 2012). For Douglas-fir, modulus of elasticity in 
turn was most affected by timber density, less by knots or 
microfibril angle (Vikram et al. 2011). However, knottiness of 
sawn timber is a major criterion for timber grade assignment 
as it influences strength properties and size (including length) 
effects (Barrett and Fewell 1990; Denzler 2007; Rais et al. 
2010; Stapel and van de Kuilen 2013). Additionally, national 
standards differ in the determination of static modulus of 
elasticity (random vs. systematic orientation) and complicate 
these dependencies. 

Previous analyses of Douglas-fir often focused on tree or 
log level with regard to wood quality and did not continue as 
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far as sawn timber quality and immediate quality parameters 
such as strength, stiffness, and knottiness of individual 
(Mitchell 1975; Carter et al. 1986; Todaro and Macchioni 
2010; Kladtke et al. 2012). A general link between growth 
and wood properties of Norway spruce was made by Houllier 
et al. (1995). This knowledge was incorporated into stand 
growth models, economic analyses, and management tools 
to investigate the impact of silvicultural practices on standing 
yield and financial return (Mitchell 1988). In this empirical 
study, we want to make precise predictions on timber quality 
disregarding volume or economic aspects. Douglas-fir 
(Pseudotsuga menziesii (Mirb.) Franco) was selected as a 
model species, but the results are of general relevance for 
softwood management. For analysis, unique Douglas-fir plots 
of the long-term forestry research plots in Bavaria, Southern 
Germany were available which were established in the 1970s 
in order to comprehend the dependencies of growth on differ- 
ent spacing. The even-aged trees represent typical commercial 
thinning trees of a wide range of managed forests. The objec- 
tives were to (1) investigate the impact of space availability on 
mechanical properties of sawn timber and to (2) quantify the 
effect of space availability on sawn timber yield by strength 
grading. The quality of side and center boards was analyzed as 
a function of initial plant density and log position (axial and 
radial) in the stem. 


2 Material and method 
2.1 Material 


A sample of 164 Douglas-fir trees (P. menziesii (Mirb.) Fran- 
co) was used in this investigation. The trees originated from 
two experimental site plots in Southern Germany which are 
part of the Bavarian network of long-term forestry research 
plots. One site is located in the north of Bavaria next to the city 
of Aschaffenburg in a low mountain range named Spessart. 
The other site is located about 15 km south of the city of 
Ansbach. Both experimental sites have been founded in the 
beginning of the 1970s. Hence, the age of all trees was about 
40 years. Within the sample, all trees were typical thinning 
trees; future crop trees were not included. At both sites, the 
experimental area is divided into smaller single plots. These 
are square shaped with an edge length of 30 m and a buffer 
zone of 5 m around the plot. Initial plant density differs within 
the single plots—1,000, 2,000, and 4,000 trees per hectare. 
Thinning was performed in accordance to a standardized 
experimental program on the basis of a guide curve described 
by Kladtke et al. (2012). Thinning type is defined as a final 
crop—tree system; thinning frequency is determined by the 
dominant height. When dominant height reached 12 m, 150 
future crop trees were selected. The guide curve gives the 
number of trees per hectare depending on the dominant height. 


High plant density improves timber strength 


In the past, thinning followed this reference curve in the 
autumns/winters of 1989, 1994, 1999, 2004, and 2010. The 
sampled trees for this research were harvested for measure- 
ment in the winter of 2010. The trees—selected by foresters— 
were typical for thinning trees according to log quality and log 
dimensions. 

Table | gives an overview of the samples. Logs of 4.1 m 
have been sawn to boards when the top diameter was at least 
24 cm without bark. Logs have been sawn using a band saw, 
with sawing pattern shown in Fig. 1. It was tried to get as many 
center boards as possible with the preferred cross-section of 
50 mm x 150 mm (n=998) as this is a standard cross-section for 
both joists loaded in bending and for lamellas used in glued 
laminated timber production. Moreover, from the outer parts of 
the stems, a number of 936 side boards with cross-section of 
30 mmx 100 mm and 30 mm 140 mm could be obtained. 


2.2 Method 


Strength grading of sawn timber can be done visually as well 
as by machines. For visual grading, the German standard DIN 
4074 (DIN 2012) is one of the most relevant in Central Europe 
(Stapel and Van de Kuilen 2013). The grading rules from this 
standard are for instance applied in countries like the Czech 
Republic, Austria, Italy, and Germany. The standard distin- 
guishes 11 properties (knots, fissures, wane, rate of growth, 
slope of grain, pith, warp, and discoloration) for which classes 
are defined. The most important criterion for the mechanical 
properties of joists and edgewise-loaded boards is knottiness 
whereby the grade is normally determined by the single-knot 
value DEK. The knot value DEK is the major ratio of the 
minimal surface diameter of an oval knot to the side width of 
the board on which the knot is visible (DIN 4074, DIN 2012). 
Only the knot value DEK was taken into consideration when 


Table 1 Sampled data, mean 
values, and standard deviation in 
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the boards were visually graded. Grades S13, S10, S7, and 
reject were evaluated. In EN 1912 (CEN 2012), grade S10 has 
been allocated into strength class C24 (Douglas-fir from Ger- 
many) according to EN 338 (CEN 2009a). Machine strength 
grading has been performed using GoldenEye-706 giving the 
dynamic modulus of elasticity (MOEdyn), a machine knot 
value, and a board density without knots calculated on X-ray 
picture of the board. 

The bending strength of the boards was predicted using the 
machine data regarding knots, density, and MOEdyn. The boards 
were assigned to European strength classes C18, C24, C30, and 
C40 (EN 338, CEN 2009a). Every single-strength class was 
graded individually using the complete data (C18/reject, C24/ 
reject, C30/reject, and C40/reject); no strength class combination 
was applied (ITT/100/10/08, CEN 2013). Strength class C40 
was graded in order to find a possible relationship between initial 
plant density (PD) and high-quality timber. The other three 
classes were graded as they are the common strength classes 
used nowadays (EN 14081-4, CEN 2009b). 

A number of 701 center boards were randomly selected and 
tested in tension according to EN 408 (CEN 2009c). All test 
pieces were of sufficient length (4.1 m) to provide a test length 
clear of the testing machine grips of nine times the larger 
cross-sectional dimension. From a preliminary test series, we 
found a constant ratio between tensile and bending strength of 
0.71 for all iPDs, thus covering the whole sawfalling range of 
wood qualities. Therefore, it is possible to use tensile strength 
as a measure for comparison also with literature presenting 
bending strength data. 


2.3 Statistical analysis 


Sawn timber quality y was described by the following vari- 
ables: visual knottiness, machine knottiness, density, 


Initial plant density iPD (trees ha ') 


brackets 
1,000 (SD) 2,000 (SD) 4,000 (SD) 
Plots Number of single plots 1] 14 8 
Trees Number 64 68 32 
DBH 330 (68) 316 (56) 305 (35) 
Height 255 (42) 270 (24) 268 (22) 
Height to DBH 79 (15) 88 (11) 88 (9) 
Logs Number 135 151 68 
Middle diameter 285 (57) 274 (46) 262 (30) 
Center boards Number 409 423 166 
Cross-section 50x 150 
Board length 4,089 (149) 4,102 (24) 4,106 (11) 
Side boards Number 407 377 152 
Cross-section mm? 30x 100, 30x 140 
Board length mm 4,071 (174) 4,091 (49) 4,100 (47) 
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Fig. 1 Ongin and location of boards in longitudinal direction (a) and 
radial direction (b); center boards were cut out from the inner part, side 
boards from the outer part of the logs 


MOEdyn, estimated bending strength, and measured tensile 
strength. An analysis was performed to verify whether sawn 
timber quality y is related to any of the following explanatory 
variables: iPD, radial position (7p), and longitudinal position 
(/p). 

The model to be fitted was then formulated as follows: 


Ving = 4 + airp + a2lp + a3iPD + 6; + by + Bie + Bim + € 


where i was the site index; j indexed the plot in site i; k indexed 
the tree on plot 7 in site 1; / indexed the log of tree k on plot 7 in 
site 7; and ao, a1, a, and a3 were the fixed effect coefficients to 
be estimated. The explanatory variables iPD (levels—1,000, 
2,000, and 4,000 trees per hectare), /p (levels—1 for butt log, 
2 for second log, 3 for third log, and 4 for top log), and rp 
(levels—center and side boards) were of categorical type and 
were dummy-coded; the reference was a main board from the 
butt log of a tree coming from a plot with an iPD of 1,000 trees 
per hectare. 

Aside from these fixed effects, a set of random effects was 
considered to match the error structure of the data. Parameters 
bi, biz, Dijx, aNd Hi, Were nested random effects at site, plot-in- 
site, tree-in-plot-in-site, and log-in-tree-in-plot-in-site levels 
(Pinheiro and Bates 2000; Zuur et al. 2009). The random 
effects were assumed to be normally distributed with zero 
mean and constant variance. These nesting level-specific var- 
lances were estimated during the model fitting process. The 
symbol ¢ represented the independent and identically distrib- 
uted random errors. 

The Akaike information criterion (AIC) was used as fitting 
criterion for selecting the best multiple linear model (Burnham 
and Anderson 2002). By means of AIC, a manual model 
simplification procedure was carried out to obtain the 
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adequate model and to remove nonsignificant parameters. 
The AIC values were calculated, and the model with the 
lowest AIC was selected and used for further analysis. It 
was tested whether this model performed better than a model 
including only random effects. 

All data were analyzed using R (R Development Core 
Team 2012) and the R package /me4 (Bates and Maechler 
2009) and languageR (Baayen 2008). The default behavior of 
[mer to optimize the restricted maximum likelihood criterion 
(REML) was used (Pinheiro and Bates 2000). Normality and 
homogeneity have been checked by visual inspections of plots 
of residuals against fitted values. 


3 Results 
3.1 Overview of data 


There was a large degree of variation in all physical and 
mechanical sawn timber properties which required power- 
ful statistical analyses. Timber quality varied between iPD. 
Table 2 illustrates the differences for relevant parameters 
of sawn timber quality. Trees growing up at high plant 
density were characterized by better timber properties. The 
knottiness of boards was lowest for iPD of 4,000 trees per 
hectare. Density, dynamic MOEdyn, and strength of 
boards increased with increasing iPD. For comparison 
purposes, some timber quality properties are mentioned 
from the literature in order to give a framework of refer- 
ence. These values were determined by small, defect-free 
specimens and cannot be directly compared to our sawn 
timber values. In average, wood density of Douglas-fir 
was given as 0.5 g/cm’, static bending modulus of elas- 
ticity as 12,500 N/mm”, bending strength as 85 N/mm’, 
and tensile strength as 95 N/mm? (Gohre 1958; Knigge 
1958; Forest Products Laboratory 1999). 


3.2 Statistical dependencies between timber properties 
and spacing as well as stem position 


As described within Section 2, all possible combinations of 
the explanatory variables rp, /p, and iPD were scrutinized with 
respect to their effect on sawn timber properties. Models were 
preferred which provided the lowest AICs. Generally, the 
models used in our analysis performed better than the accord- 
ing models that included only random effects. Timber quality 
describing parameters—knottiness (either visual or by ma- 
chine), density, MOEdyn, estimated strength, and tensile 
strength—were best explained by all three predictors (rp, Ip, 
and iPD). Although the trees originated from two different 
plots in Southern Germany, site was not of much relevance. 
The rp has significant influence if the center and side 
boards were compared: boards sawn from the outer part of 


High plant density improves timber strength 


Table 2 Sawn timber properties depending on iPD, mean value, standard deviation in brackets 


Center board Knottiness DIN 4074, DEK* 
Machine 
Density 
MOEdyn 
Estimated bending strength 
Tensile strength* 
Side board Knottiness DIN 4074, DEK 
Machine 
Density 
MOEdyn 


Estimated bending strength 
Tensile strength 


621 
Initial plant density iPD (trees ha ’) 
1,000 (SD) 2,000 (SD) 4,000 (SD) 
- 0.53 (0.18) 0.47 (0.15) 0.42 (0.16) 
— 2,639 (497) 2,312 (519) 2,126 (491) 
kg m° 461 (42) 466 (43) 474 (42) 
N mm” 10,509 (1,898) 11,805 (2,291) 12,578 (2,232) 
N mm” 24.3 (10.3) 315 (11.5) 35.8 (10.6) 
N mm~ 17.3 (6.9) 22.5 (9.2) 25.7.7) 
- 2,421 (616) 1,979 (631) 1,684 (568) 
kg m° 494 (44) 504 (44) 511 (47) 
N mm” 11,622 (2,171) 13,443 (2,640) 14,327 (2,605) 
N mm” 32.4 (10.5) 41.5 (11.7) 46.1 (11.0) 
N mm ” — — — 


“Knottiness according DIN 4074 (DEK) and tensile strength were only determined on a smaller sample size: 92, 311, and 98 in increasing order of iPD 


the stem were characterized by higher quality (see Table 3); 
knottiness was lower, density (35 kg/m°), MOEdyn (1,513 N/ 
mm”), and estimated bending strength (9.6 N/mm”) were 
higher. 

Along the height of the tree (Jp), Douglas-fir showed de- 
creasing sawn timber quality. Sawn timber from a height of 
about 11 m above ground (/p=3) had a significantly lower 
(—21 kg/m°, Table 3) density (444 kg/m°) than sawn timber 
from the butt log (465 kg/m”). There was also a difference in 
MOEdyn and strength of sawn timber between the different 
longitudinal positions in the stem. Sawn timber MOEdyn 
from a butt log (Iy=1) was 507 N/mm” higher than from the 
fourth log (Jp=4); estimated bending strength was 5.1 N/mm? 
higher. The measured tensile strength of jp=3 was 2.5 N/mm? 
lower than that of Jp=1 (see Table 3). 

High iPD improved timber strength properties positively 
(see Table 3). Except for wood density, all other properties 
were significantly different between iPD=1,000 and iPD= 
2,000 trees per hectare as well as between iPD=1,000 and 
iPD=4,000 trees per hectare. The differences between iPD= 
2,000 and iPD=4,000 trees per hectare were not significant at 
a probability level of 0.05. Knottiness—visual and machine— 
significantly decreased with iPD for both kinds of boards. The 
parameter MOEdyn—as one of the best predictors for 
strength—was positively affected by iPD. For estimated and 
measured strength, similar results were observed. For in- 
stance, tensile strength values increased by about 6 N/mm” 
when iPD was doubled from 1,000 to 2,000 trees per hectare. 
There was only a slight, but not significant, difference of 
tensile strength between trees growing in density of iPD= 
2,000 and iPD=4,000 trees per hectare. Table 4 shows the 
yield for different grades depending on the iPD. Only boards 
of the butt logs were considered. The kind of grading— 


whether visual or machine grading—did not influence the 
general trend: Yields were the highest for boards cut from 
trees with high iPD for all strength classes. For instance, none 
of the center boards from trees growing up at iPD= 1,000 
reached C40, and only 10 % of the center boards of the lowest 
iPD class met the requirements of C30; yield rose to threefold, 
1.e., 32 % for iPD=4,000. The yield of machine-graded center 
boards in C24 rose by about 40 % from the lowest (50 %) to 
the highest (89 %) iPD. Side board yield also achieved higher 
levels for each grade. It increased from 84 to 99 % for side 
boards of grade C24. The yield of visually graded timber from 
iPD=4,000 was more than twice as high as from iPD=1,000. 
In comparison to machine-graded timber, visually graded 
timber reached only half the yield. For Douglas-fir, visual 
grading seems to be not efficient (Roblot et al. 2008). 


4 Discussion 
4.1 High variation of timber quality of Douglas-fir 


The sampled wood represented young thinning trees of log 
dimensions which were usually sawn to structural timber. The 
low quality was partly caused by the age of 40 years and the 
high proportion of juvenile timber; for Douglas-fir, about the 
20th ring approximates the point of transition from juvenile to 
mature wood based on the pattern of pith-to-bark wood den- 
sity (DiLucca 1989; Barrett and Kellogg 1991; Kennedy 
1995). Juvenile wood is characterized by a lower percentage 
of late wood, lower density, shorter tracheids with larger 
microfibril angle and grain angle, and consequently lower 
strength and stiffness properties (Bendtsen 1978). Glos et al. 
(2006) used a very similar material in a study of Douglas-fir in 
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Table 3 Effect of radial position (rp), longitudinal position (/p), and iPD on sawn timber properties regarded as knottiness, density, MOEdyn, and 


strength 
Knottiness 
Visual Machine 
ag 0.51 2,501 
0.03 61 
a), Vp Center —a 
Board = 
Side —308 b 
Board at 
a>, lp 1 —a —a 
2 +0.003 a +146 b 
0.014 oy 
3 +0.069 b +286 c 
0.019 39 
4 +0.036 ab +391 ¢ 
0.030 72 
a3, iPD 1,000 —a —a 
2,000 —0.06 b —394 b 
0.02 76 
4,000 —0.10b 507 b 
0.03 89 
b; o 0.76x 10° 10510 
bi o 1.33x10° 14.5x10° 
dix o 0.98x10° 61.8x10° 
Digs a 1.84x 10° 10.8 10° 
E om MI 202.5 x 10° 
n 701 1,934 


Density MOEdyn Estimated bending Tensile 
strength strength 
kgm° N mm” N mm” N mm ~ 
465 10,588 25.5 19.1 
7 323 1.3 i 
=a —a =a 
+35 b +1,513 b +9.6 b 
68 0.4 
=a = -d —a =d 
-14b —90 ac —0.9 b -1.4b 
96 0.5 0.6 
—21¢ 509 b -3.4¢ 2.5 b 
142 0.7 0.8 
—24¢ —507 bc =—§,1¢ —1.9 ab 
252 es 13 
=a —a —a —a 
+19 b +1,909 b +92 b +6.2 b 
410 1.7 1.4 
+13 ab +1,944 b +10.4 b +7.1b 
10 486 2.0 17 
15 0.3x10° 0.0 0.0 
186 3.8x10° 5.9 3.7 
1,031 25.4x10° 43.8 23.5 
52 2.0x 10° 32 4.0 
517 2310 64.6 36.6 
1,934 1,934 1,934 701 


Within a category rp, /p, and iPD, numbers within the same column are significantly different at the 0.05 probability level if they are different in all letters 


(a—c); standard errors in italics 


Germany, and found a mean tensile strength of 24.5 N/mm” 
and a mean MOEdyn of 12,400 N/mm? (n=315). Sauter 
(1992), on the other hand, used for his investigation about 
100-year-old trees from Germany with higher sawn timber 
quality and took a sample of 250 boards that had a mean 
MOEdyn of about 15,000 N/mm? (our sample—12,100 N/ 
mm) and a mean tensile strength of over 42 N/mm” (our 
sample—20.8 N/mm”, only main boards). The stiffness and 
strength properties of this study apparently represented the 
lower tail of the distribution of mechanical properties. 


4.2 Influence of radial and longitudinal board position 
on wood quality 


The results showed a consistent effect of the distance to pith 
on all mechanical properties studied. Side boards were of 
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higher quality than center boards. The estimated bending 
strength of side boards was almost 10 N/mm” higher than that 
of center boards and hence immediately affected sawn timber 
yield. Our result went in line with Glos et al. (2006) who 
separated their sample into boards with and without pith: 
Tensile strength of boards with pith was 21.5 N/mm? as 
compared to 27.6 N/mm”, and MOEdyn was 11,500 N/mm? 
as compared to 13,400 N/mm”. Previous investigations con- 
firmed the change of timber characteristics in radial direction: 
for sawn timber of Norway spruce (Brichert et al. 2000; 
Oscarsson et al. 2011) and for small, defect-free specimens 
of young Douglas-fir (Langum et al. 2009). 

In general, sawn timber quality of Douglas-fir decreased 
with increasing distance to ground (Knigge 1958; Barrett and 
Kellogg 1991). We observed a gradual decline of MOEdyn 
with increasing longitudinal position that was also reported by 


High plant density improves timber strength 


Table 4 Sawn timber yield (in percent) of boards depending on iPD, only 
boards from butt logs 


Grading  Strengthclass Initial plant density iPD 
(trees ha ') 
1,000 2,000 4,000 
Yield (%) Center board 
Visual S10+ (C24) 26 os 54 
Machine C18 83 23 99 
C24 50 7S 89 
C30 10 phe 32 
C40 0 2 - 
Side board 
Machine C18 96 98 100 
C24 84 95 99 
C30 33 66 TH 
C40 1 7 14 


Langum et al. (2009). Similar height-dependent MOEdyn 
patterns were also observed for Scots pine (Lindstrom et al. 
2009) and for Norway spruce (Bruchert et al. 2000). However, 
in our analysis, the difference was low between boards cut out 
from butt and second log. The significant difference was about 
500 N/mm” between MOEdyn of the butt log and the third 
log. Sawn timber density got lower with increasing height 
position and confirmed results of Barrett and Kellogg (1991) 
who studied second-growth Douglas-fir 2 x 4 lumber from 46 
to 60 years. Strength behaved in accordance to Langum et al. 
(2009)—decreasing strength values with increasing vertical 
position within a tree—although they investigated clear (de- 
fect free) wood specimens. 


4.3 Influence of initial plant density on wood quality 


Initial plant density iPD had significant influence on sawn 
timber quality for structural use. Within an iPD range of 
1,000, 2,000, and 4,000 trees per hectare, it clearly determined 
knottiness, density, MOEdyn, and strength. For all strength 
classes under consideration, the highest plant density achieved 
the highest yield. Previous studies assessed wood quality at 
the stage of logs and showed that iPD was an important factor 
which affected the size of branches (Brazier and Mobbs 
1993), specifically in the butt and second log of Douglas-fir 
(Robbins 2000). Pfister et al. (2007) also found a positive 
correlation between the largest branch diameter of log and 
initial density in a 33-year-old Norway spruce plantation in 
southern Sweden. The added value of this study was that it 
clearly demonstrated the dependence of sawn timber quality 
on iPD both for center and side boards: It confirmed that 
knottiness decreases with increasing iPD and at the same time 
proves that MOEdyn and strength both increase. Some 
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research has also been conducted in previous work on stiffness 
of structural timber, but quantitative data is limited: Based on 
Norway spruce in Sweden, Johansson (1992) concluded that 
the quality of machine stress-graded timber (deflection) in- 
creased in response to a decrease in spacing. For Sitka spruce, 
increased planting distances influenced stiffness properties of 
sawn timber (measured by a deflection grading machine) 
negatively (Brazier and Mobbs 1993). In our study, also wood 
density increased with iPD, albeit not consistently and only 
from 1,000 to 2,000 trees per hectare at a significance level of 
0.05. This was in accordance to the findings of Johansson 
(1993) who reported a significant but only small effect of little 
practical importance for a Norway spruce stand of similar age 
and structure. Fujimoto and Koga (2009) reported also on an 
unclear dependency between initial spacing and wood density 
of 22-year-old larch, but they observed that closer spacing 
results in a faster increase in wood density from the pith 
outward than for more widely spaced trees. They concluded 
that initial plant density directly influences the age of transi- 
tion from juvenile to mature wood. Low planting distances 
had also a negative influence on the proportion of juvenile 
wood of Sitka spruce (Brazier and Mobbs 1993). 

In terms of quality, the ideal silviculture suggests high 
stand density management. Nevertheless, foresters have not 
reacted noticeably to these findings because of higher costs 
involved in increasing initial stocking, the extra thinning 
which may be required, and the additional rotation age. The 
results of this study give foresters the option not to focus on 
the optimization of volume alone (current practice) but also on 
that of wood quality. In practice, however, foresters barely 
participate on the added end-product value. As foresters are 
paid by product volume and not by sawn wood quality, there is 
no direct need for foresters to improve their operations such 
that higher wood qualities are obtained. Logs grading is also 
promising, but logs are valued mainly on the basis of volume 
and not on the basis of the wood quality sawn from them. 
Sawn timber of high-strength classes achieve higher prices on 
the market, and sawmills are able to estimate the end-product 
quality as well as the value very accurately by machine 
strength grading. This knowledge could be considered in the 
current roundwood pricing system to give incentives for a 
more quality-oriented forest management. 


4.4 Limitations of the transferability to other softwood species 
and to older trees 


Wood density of other conifer species such as spruce or pine is 
likely to depend on iPD in the same manner as that of Douglas- 
fir: for all conifers, high iPDs result in low radial growth rates, 
the proportion of latewood within annual rings increases, and 
finally wood density increases (Kollmann 1951), but probably 
the presence of juvenile wood obscures this rule. The depen- 
dencies between iPD and the other sawn timber properties can 
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be transferred to both conifers and deciduous trees, as growing 
patterns are similar with respect to space availability: The iPD 
influences the branchiness of trees and consequently the knot- 
tiness of boards that affects mechanical properties directly, 
because wood fibers are deviated by local defects as knots. 
Larger knots cause higher timber grain deviation, a lower 
stiffness and especially lower strength. For Douglas-fir, wood 
density seems to explain timber properties as stiffness and 
strength more accurately than microfibril angle (Lachenbruch 
et al. 2010; Vikram et al. 2011) or knottiness (Vikram et al. 
2011). Even though some literature emphasized that density is 
more important than knottiness with respect to mechanical 
properties, our research indicated that this is only valid for 
modulus of elasticity, but clearly not important for tensile 
strength of structural size of Douglas-fir. Our investigation 
confirmed that MOEdyn showed a stronger dependency on 
density (7°=0.52) than on machine knottiness (7 =0.37) for all 
1,934 boards. The 7* between tensile strength and density was 
0.18 and between tensile strength and knottiness was 0.39 (n= 
701). Timber quality of larger (older) Douglas-fir trees is likely 
to depend on iPD as competition prevails and is likely to exert 
continuous and positive influence on the yield of high-quality 
boards due to wood density as well as thickness of branches 
and knots. The dependence of wood quality on iPD is evident 
for the lower part of the stem in particular. Hence, it may be 
expected to even persist far beyond the juvenile phase of stand 
development to the stage of final cutting. 


5 Conclusion 


The analysis aimed to contribute to the facilitation of a pro- 
duction chain that ranges from the management of softwood 
forests to the decision making of sawmills. It clearly demon- 
strated that foresters can influence end-product quality by 
controlling initial space availability of trees. The sample rep- 
resented young thinning trees of an age of 40 years whose logs 
were mostly processed to sawn timber. The study showed that 
increased radial and decreased axial position of the boards 
strongly determined mechanical sawn timber properties posi- 
tively. A high initial number of trees per hectare resulted in 
lower knottiness, higher timber density, and higher stiffness 
and strength. The influence on boards’ knottiness, MOEdyn, 
and strength was observed in particular between an initial 
plant density of 1,000 and 2,000 as well as 1,000 and 4,000 
trees per hectare. 

Our findings provide additional facts that help to derive the 
appropriate stand density in order to put into effect a specific 
management objective. Many long-term spacing and thinning 
experiments showed that a stand density reduction below 20- 
30 % of the maximum stand density can severely reduce the 
mean annual increment and total volume yield of forest stands 
(Pretzsch 2005). Wide spacing and early thinning, on the other 
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hand, can significantly accelerate diameter growth and shorten 
the rotation period. The findings that wood quality reductions 
are associated with low densities provide new arguments for 
keeping higher density levels for achieving both high-quality 
and quantity of wood production. Future forest management 
plans should provide a balanced mixture of strategies 
concerning wood quality and volume as well as risk manage- 
ment on the level of stands. 

Our results call into question a low-density management on 
fertile sites in particular where growth rates are high a priori, 
and strong thinning results in even higher growth rates, wider 
tree rings, lower wood density, and severe yield reduction. On 
poor sites, in contrast, the effect of thinning on quantity and 
quality of wood production might be less disadvantageous. 
Future evaluations can be based on the shown trade-off be- 
tween density and wood quality in order to weigh and opti- 
mize the pros and cons of density management by integration 
into exiting growth models and scenario analyses. 

Information about treatment effects on quality and 
quantity of the produced timber in a given stand develop- 
ment phase or even over the whole rotation were so far 
nearly missing. Thus, an overarching optimization of eco- 
logical, production, and economical results of Douglas-fir 
management in forest practice was not possible. This will 
be remedied when the revealed relationship between stand 
density management and sawn timber properties of Doug- 
las-fir will be integrated in individual tree growth models. 
This enables scenario analyses about the effect of stand 
density on a broad range of output variables including 
wood quality and quantity. Here, we provide the indis- 
pensable knowledge for such a model completion. In 
contrast to the quickly changing technical and economical 
frame conditions, the revealed between-tree and within-tree 
wood allocation patterns are time invariant and indispens- 
able for further evaluations of treatment effects at tree and 
stand level. To our state of knowledge, the analyzed 
Douglas-fir stands are the only available samples of this 
species in Europe where both timber quantity and quality 
have been recorded since stand establishment and provide 
a unique database for model development and calibration. 
Finally, further investigations should focus on older 
(future) trees or model the quality development with in- 
creasing tree age. However, due to the reduced knottiness 
and its positive influence of the yield of side boards, it is 
likely that high initial plant density also promotes timber 
quality of larger trees. 
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